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Abstract 

The radiative decays of the phi meson are known to be a good source of information 
about the ao(980) and /o(980) scalar mesons. We discuss these decays starting from 
a non-linear model Lagrangian which maintains the (broken) chiral symmetry for the 
pseudoscalar (P), scalar (S) and vector (V) nonets involved. The characteristic feature 
is derivative coupling for the SPP interaction. In an initial approximation which models 
all the scalar nonet radiative processes together with the help of a point like vertex, 
it is noted that the derivative coupling prevents the ao and fo resonance peaks from 
getting washed out (by falling phase space). However, the shapes of the invariant 
two final PP mass distributions do not agree well with the experimental ones. For 
improving the situation we verify that inclusion of the charged K meson loop diagrams 
in the model does reproduce the experimental spectrum shapes in the resonance region. 
The derivative coupling introduces quadratic as well as logarithmic divergences in this 
calculation. Using dimensional regularization we show in detail that these divergences 
actually cancel out among the four diagrams, as expected from gauge invariance. We 
point out the features which are expected to be important for further work on this 
model and for learning more about the puzzling scalar mesons. 



2 



1 Introduction 

Recently, there have been a number of important experimental studies |T] of the rare radiative 
decays of the 0(1020) vector meson: — > tttcj and — > 71777. These decays seem to be 
dominated by the production (and subsequent decay) of the scalar mesons, a (980) and 
/o(980) according to — > fo, a + 7 and hence are generally considered to provide valuable 
information about the puzzling light scalar mesons|2j of low energy QCD. 

The theoretical analysis of this type of decay was initiated by Achasov and Ivanchenko 
|3| and followed up by many others [3]. The starting point was the observation that the 
meson decays about 50 per cent of the time into K + K~ . Since this final state can easily 
annihilate to produce either an fo or ao together with an emitted photon, it is rather natural 
to consider charged i^-meson loop diagrams to describe the process. Similarly the meson 
is observed to decay about 15 per cent of the time to np or 7r + 7r~7r° so one expects some non 
resonant background which is likely to include the emission of a pion with a virtual p which 
subsequently decays into ify (and similar diagrams leading to a ir rj^ final state). 

The varied calculations along these lines lead to results which more or less agree with 
experiment. Of course it is desirable to fine tune this agreement, both to reflect the expected 
improved accuracy of new experiments as well as to improve our understanding of strong 
interaction calculations. Here we will focus on some technical points, which do not much 
change the previous results but may be of interest for future more ambitious calculations as 
more experimental data become available. Mainly, we will require that the amplitudes all be 
computed from a chiral invariant Lagrangian (containing usual quark mass induced breaking 
terms). This is a symmetry of nature apparently so it is desirable to calculate in this way 
even though the spontaneous breakdown of chiral symmetry (in the absence of quark mass 
terms) means that, especially away from thresholds, one can often get reasonable predictions 
by not explicitly taking it into account. 

Two approaches are commonly employed to implement the chiral symmetry in the ef- 
fective Lagrangian framework. In the linear sigma model approach, scalar partners of the 
pseudoscalars are introduced. In the non-linear sigma model approach, one initially deals 
with pseudoscalars only, the scalars having been essentially "integrated out" . The character- 
istic feature of the non linear model is the appearance of derivative type interaction terms as 
opposed to non derivative type interaction terms in the linear model. Nevertheless, the non 
linear model is often more convenient to use. For example, the celebrated result |3] for near 
threshold pi pi scattering arises in the linear model from a delicate cancellation of two rather 
large terms. On the other hand it arises directly from a simple single term of the correct 
characteristic strength in the non linear model. In the present paper we shall deal with the 
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non linear model approach. Since vector and scalar mesons are also involved in the processes 
of interest we will add these to the non linear Lagrangian of pseudoscalars in a conventional 
way. Such a formulation essentially implements vector meson dominance automatically for 
processes involving photons. 

We shall restrict our attention further here to processes of the type — > 7 + virtual scalar 
where the virtual scalar (either or / ) subsequently decays to two pseudoscalars. First 
we shall consider a possible non - K + loop contribution to this process. We previously jH] 
studied this by introducing an effective strong VVS (vector-vector scalar) interaction based 
on an analogy to the effective VVP (vector vector pseudoscalar) interaction used many years 
ago [7j to study analogous processes like a; (782) — ► tt°j. This might open the possibility of 
understanding properties of the whole nonet of scalars at once. Especially, it might shed 
some light on the composition of the light scalar nonet; whether the light scalar mesons 
are composed of one quark and one anti-quark (2-quark picture) or two quarks and two 
anti-quarks (4-quark picture). 

In the present paper, we point out an interesting effect. If a non derivative SPP type 
interaction were to be used there would be a strong tendency for the decreasing phase space 
to wash out the predicted scalar meson peak in, for example, dT(<j) — > TT°i]'~f) / 'dq. Here q 2 is 
the invariant squared mass for the 7ir) system. On the other hand, the use of a derivative type 
SPP interaction, as is required for chiral symmetry in the non linear sigma model approach, 
restores the peak. (There is not necessarily any contradiction with the expectation that the 
same physics near threshold should be expressed by suitably generalized linear and non linear 
models. One expects the linear model description to include additional terms). However, 
we notice that there is experimentally more enhancement of the scalar peak than can be 
accounted for by this mechanism. Thus we are led to also consider the usual K + loop 
diagram in our approach. 

As mentioned, the K + loop diagram has been considered by many authors OH]. We 
can not basically change the well established results. However we note that the effect of the 
derivative couplings will also sharpen the scalar peak. Actually, the derivative SPP couplings 
result in quadratic as well as logarithmic divergences and an additional diagram. It has been 
found jS] that such "unpleasant details" of the calculation can be circumvented by assuming 
gauge invariance. Specifically, gauge invariance requires that the amplitude for — > photon + 
scalar be proportional to e^ejf (d^p-k— p^k^), where e v and p are respectively the polarization 
and momentum four vectors of the (f) meson while e and k correspond to the photon. Then 
it is only necessary to calculate the coefficient of the p^ky term, which eliminates the need to 
calculate two diagrams and worry about divergences actually cancelling each other. Of course 
it would be nice to regulate all the diagrams and verify in detail how the cancellations take 
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place. We have carried out this somewhat lengthy task using the dimensional regularization 
scheme and will give details in the present paper. 

In section |2J we first present the chiral Lagrangian of pseudoscalars, vectors and scalars 
which will be used for the subsequent calculations. Our initial motivation, described in 
Ref. 0, was to relate all the decays of the types S— > 77, V — >Sj and S — >V7 to each other 
by using a simple effective point like interaction. We next consider the 0(1020) decays into 
7r°rj and 7r°7r° proceeding respectively from intermediate ao(980) and /o(980) resonances in 
this simple model. It can be seen that the spectrum shapes for large q are not as sharply 
peaked as the experimental data indicate. 

In section |3J we calculate the form of the charged K meson loop contributions to these 
two decays using a non-linear chiral Lagrangian which maintains the chiral invariance when 
vectors and scalars as well as pseudoscalars are included. The extension to include photon 
interactions is given. It is noted that individual diagrams contain quadratic as well as 
logarithmic divergences. A careful treatment using the dimensional regularization scheme 
shows that these divergences both cancel leaving a finite answer. 

In section |Uwe study the spectrum shape of the i^-loop contributions to these decays. We 
find that this has a characteristic shape which does in fact agree with experiment, suggesting 
that the dynamics of the K loop plays an important role. 

Section El contains a brief summary. Some discussion will be given on the status of the 
present program and related future work. 

2 VVS type contributions to <\> — > 7T°?77 and — > 7r°7r°7 

Our calculation is based on a standard non-linear chiral Lagrangian containing, in addition 
to the pseudoscalar nonet matrix field 0, the vector meson nonet matrix p M and a scalar 
nonet matrix field denoted by N. Under chiral unitary transformations of the three light 
quarks; ^l.r — > C/l,r • Ql,k, the chiral matrix U = exp(2z0 / F n ) , where F n ~ 0.131 GeV, 
transforms as U — > ■ U • U^. The convenient matrix K(Ui,,Ur, 0) is defined by the 
following transformation property of £ (U = £ 2 ): £ — ► ■ £ • = K • £ • U^, and specifies 
the transformations of "constituent-type" objects. The fields we need transform as 



N -»• K ■ N ■ K ] , 

p, - K ■ p, ■ + 4^ • , 



FtMu(p) = d^py - d v p lx - ig [p^ , p v ] -> K ■ ■ K ] , 
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where the coupling constant g is about 4.04. One may refer to Ref. ^U] for our treatment 
of the pseudoscalar-vector Lagrangian and to Ref. jTTj for the scalar addition. The entire 
Lagrangian is chiral invariant (modulo the quark mass term induced symmetry breaking 
pieces) and, when electromagnetism is added, gauge invariant. The U(3)l x U(3)r invariant 
portion of the effective Lagrangian reads: * x 

Co = - ^Tr(^) - ~Tr(F M „(p)F Mi ,(p)) 

- l -Tv{V,NV,N) - ||tt [(gp, - v,f 

- aTr(NN) - C Tr(iV)Tr(iV) 

+ F*[Ae abc e def N*(p,y b (p,)l + 5Tr(iV)Tr( W ) 

+ CTr(Np li )Tr(p ll ) + DTr(AT)Tr(^)Tr(p M )] (2.2) 

where V^N = d^N — igp^N + iNgp^.* 2 Furthermore v^,p^ = (i/2)(£d fl ^ ± £*<9 M £), where 
£ = U 1 / 2 . These terms include the parameters m^,a, c, A, B,C and D. More details about 
the evaluation of these parameters are discussed in Refs. [TO] and [TU] . 

It should be remarked that the effect of adding vectors to the chiral Lagrangian of pseu- 
doscalars only is to replace the photon coupling to the charged pseudoscalars as, 



ieAy, Tr [Q<p > 

kgF 2 Tr(Q P „) 



oAfj, 



k N 



+i\l-~\Tr[Q<!>d fl , 



+ •••, (2.3) 



where is the photon field, Q = diag(2/3, — 1/3, — 1/3) and k = {jp-j with m v ~ 
0.76 GeV. The ellipses stand for symmetry breaking corrections. We see that in this model, 

# 1 This Lagrangian can be rewritten within the framework of the hidden local symmetry (HLS) |12l I13| . 
The method of including vector mesons used in this paper based on the proposal in Refs. |14| is equivalent 
to that based on the HLS approach at tree level. JS] When we consider the vector mesons inside the loop, 
the two approaches might have some differences. In the present analysis, however, we will consider the loop 
corrections from only the kaon, which provides a large enhancement to the <f> radiative decay amplitude. All 
other loop corrections from vector mesons are naturally expected to be small. In this sense, the method used 
in this paper is completely equivalent to the recently developed method ^3] used in the HLS. Note that the 
scalar mesons have not been included inside the loop in either approach. 

# 2 One could also use for the covariant derivative, the combination cgp^ + (1 — cjv^ with c being an 
arbitrary constant. In any case, there are a few more terms such as tr((<?p M — v^) N (gp^ — u^) N) and 
tr {^{gPf! — v^) 2 N 2 ^j , which include the same number of derivatives. We note that the above extra terms as 
well as the interaction terms from the covariant derivative do not contribute in the present analysis, where 
we are considering the processes related to only one scalar meson. 
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Sakurai's vector meson dominance simply amounts to the statement that k = 2 (the 
KSRF relation [TH]). This is a reasonable numerical approximation which is essentially stable 
to the addition of symmetry breakers fUJ and we employ it here by neglecting the last 
term in Eq. (|2~3|) . 

The proposed effective SVV type terms in the effective Lagrangian are 

Csvv = (3Ae abc e a ' b ' c ' [F^(p)] a a , [F^(p)] b b , N c c , 
+ P b Tt[N}Tt[F, v (p)F imu ( P )] 
+ /3 c Tr[NF flv (p)]Tr[F ftv (p)} 

+ (3 D Tr[N}Tr[F^(p)}Tr[F fiV (p)} . (2.4) 

Chiral invariance is evident from Eq. (J2.1)) and the four flavor-invariants are needed for 
generality. (A term ~ Tr(FFN) is linearly dependent on the four shown). Actually the (3d 
term does not contribute in our model so there are only three relevant parameters (3a, (3b 
and (3q- 

2.1 ao(980) production 

The Feynman diagram for the contribution from the new VVS terms to the decay process 
<f>{p, eV ) ~^ 7r °(Qi)v('l2}l(k,€) is shown in Fig. [TJ Note that the photon is produced through 




Figure 1: Feynman diagram for cf>(p,e v ) — > ir° (qi)r)(q2)j(k, e) using an effective VVS term 
its mixing with vector mesons according to Eq. ()2.3)1 . The Feynman amplitude is 
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e(qi ■ g 2 )n ( r ) [(P • *0(e V • e) - (p • e)(* • e v )] , (2.5) 

where 

^ = ^^ (g 2 )7a^. (2.6) 

Here is given in terms of the coefficients of Eq. (J2.4)) and a scalar mixing angle in Eq. (8) 

of Ref. jH] and will be considered, for generality, a single parameter. Furthermore we use the 
simple a propagator: 

D ao (q 2 ) = — =- • (2.7) 



Also, q is the positive quantity: 



q 



(po - k ) 2 - (P - k) 



1/2 



(2.8) 



Finally, the SPP type coupling constant in Eq. (J2.6|) as well as others needed in this paper 
are defined from the Lagrangian density: 

- ?sZ£e%d ll K-dpK + - ^Mf d,K-d,K + + ■■■ . (2.9) 

The relations between these coefficients to A, B, C, D in Eq. ()2.2j) are given in Appendix C 
of Ref. [TI]. The "g-distribution" dT(<f) — > Tr°rpf)/dq is expressed as 

dT{(j) -> n°T]-f) a / M|-g 2 \ 3 / [g 2 - (m, + Try) 2 ] [g 2 - (m, - m n ) 2 ] 
dq ~ 768tt 2 ^ M / V g 2 

x (q 2 - ml - mlf \Y^\ 2 . (2.10) 

Discussion of the phase space integral is given, for example, in Ref. (2Hj. 

Now let us see how well we can fit the experimental data on the tt°i] invariant mass 
distribution in this model. We will use the inputs: m ao = 984.7 MeV (from the PDG 
table [201); r a = 70 MeV (from [2]]); J ao ir V = -6.80 GeV^ 1 (from |2HE2|). 

Let us perform two types of fits for obtaining the best value of (assuming g to be 
fixed at the value 4.04): 

(I) use the data for all values of q = m^o, , (2-H) 

(II) use the data for > 850 MeV . (2.12) 



The results are 
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(I) C;° = 3.7±Q.lGeV~ 1 , xVd.o.f = 41/(32 - 1) , 

(II) =3.6±0.1GeV- 1 , X*/d.o.f = 32/(17 - 1) . (2.13) 

Figure El shows the resulting plots of dB(<p — > ■K°r] r ))/dq together with the experimental data. 
Note that, since only the combination 7 ao7r r ? C?° jg appears in our fitting procedure, the best 
fitted curve will not change even if we allow the values of 7 ao7 r» ? and g to vary. 




750 800 850 900 950 1000 



Figure 2: dB((f> — > n°r]j)/dq x 10 7 (in units of McV 1 ) as a function in the 7r°-?7 invariant 
mass q = m^o v (in MeV). Solid line shows the ao contribution with the best fitted value 
C?° = 3.7GeV~\ and the dashed line shows that with C£° = 3.6 GeV -1 , Experimental data 
indicated by white diamonds (O) are from the SND collaboration in Ref. [231, and those by filled 
triangles and filled diamonds are shown in Ref. [21] extracted from the KLOE collaboration in 
Ref. |2HJ- 

This model gives a poor fit to the experimental data in the energy region above 950 MeV. 
One possibility is that the fit may be improved by raising the mass of ao above 984.7 MeV. 
Actually, Ref. j23] gives the best fit value as m ao = 995lfg MeV. Let us then fit the a mass 
together with value of C^°. The results are 

(I) Q° = 4.0 ± 0.1 GeV _1 , m ao = 993.2 ± 2.8 MeV , xVcLo.f = 39/(32 - 2) , 

(II) Cf = 3.9 ± 0.1 GeV" 1 , m ao = 990.4 ± 2.5 MeV , X*/d.o.f = 31/(17 - 2) . (2.14) 

Note that the best fit value of m aQ in case (II) is very close to the values shown in Ref. [23J. 
In Fig. El we plot dB(<f) — > ■K°r]^)/dq together with the experimental data. This figure shows 
that it is still difficult to reproduce the experimental data in the energy region above 950 MeV 
in the present model even if one allows the ao mass to vary. 

For comparison with the chiral symmetric case, we will now investigate the effect of 
using non-derivative coupling at the ao7r°r] interaction vertex. This amounts to multiplying 
Eq. (j2~TUj) by the factor: 
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Figure 3: dB(cj> — > n°r]^)/dq x 10 7 (in units of McV 1 ) as a function in the 7r°-?7 invariant 
mass q — m^a rj (in MeV). Solid line shows the ao contribution with the best fitted values 
Cf = 4.0 GeV" 1 and m ao = 993.2 MeV, and the dashed line shows that with = 3.9 GeV" 1 
and m ao = 990.4 MeV. Experimental data are as in Fig. |2 

(m 2 aQ -ml- m 2 ) 2 



(q 2 — m 2 — m 2 ^ 2 



(2-15) 



-it j 

which has the effect of deemphasizing the high q region. It yields 



(I) Cf =2.13 ±0.07 Ge^ 1 , x 2 /d.o.f = 113/(32 - 1) , 

(II) Cf = 2.68 ± 0.08 GeV- 1 , v 2 /d.o.f = 67.9/(17 - 1) . (2.16) 

Furthermore, Fig. 0] shows the plot of dB(<f) — > 7r°?77) / ' dq together with the experimental 
data. Comparing this figure with Fig. |2]and the results in Eq. ()2.16|) with those in Eq. (|2.13J) 
indicates that the derivative coupling model gives a better fit. The non derivative coupling 
factor clearly seems to wash out the resonance peak. 

2.2 / (980) production 

The treatment of the decay <p — > 7r°7r°7 assuming only the VVS type interaction where S is 
identified as / (980) and subsequently decays to the two neutral pions, proceeds in a similar 
manner. Again it is found that the use of a chiral symmetric derivative type interaction is 
to be preferred because it does not wash out the scalar resonance peak. However the overall 
fit to the 7T7T invariant mass distribution is not good, again suggesting that the VVS type of 
contribution is not the dominant one. In this case, dT{cf) — * TT°7r°^f)/dq is given by, 



dq 15367T 2 

where 




q 2 - Ami (q 2 - 2ml) 2 \Y^\ 2 , (2-17) 
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Figure 4: dB(cj> — > Tr°r]^)/dq x 10 7 (in units of MeV -1 ) as a function in the 7r°-77 invariant 
mass q — m,o, (in MeV). Solid line shows the ao contribution in the non-derivative coupling 
model with the best fit value = 2.13 GeV -1 and the dashed line shows that with = 
2.68 GeV -1 Experimental data are as in Fig. [5] 



Y^ = ^D f0 (q 2 )V2 lf0 ^. (2.18) 



The /o propagator is: 



and we will use the mass of the /o(980) to be 987 MeV. The coupling constant 7/ 07r 7r is 
related to the width of f as 



-1 " 'T/OTT 

fo 



lm *(ml-2ml) 2 . (2.20) 



m% 



647r m/ \ 

In Ref. [26J a treatment of tot scattering suggested r/ ~ 64.6 MeV and correspondingly 
It/ottttI ~ 2.25 GeV -1 . Considering both Tin and irK scattering, 7/ 07r7r ~ 1.47 GeV -1 and 
r/ ~ 27.6 MeV were determined in Ref. 

Using for example |7/ 07r7 r| = 1-47 GeV -1 let us next fit the value of CP to the experimental 
data. Furthermore, to avoid any possible confusion with an expected low energy contribution 
from the a we shall use experimental data only in the region 

m^o > 850 MeV . (2.21) 

This yields 

C$ « 9.3GeV -1 , X 2 /d.o.f. = 101/(17 - 1) . (2.22) 
In Fig. we show the resultant fo contribution together with the experimental data [23 1 
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Figure 5: dB((j) — > 7T 7r°7)d<7X 10 8 (in units of MeV -1 ) as a function of the dipion invariant mass 
q = m n o„o (in MeV). Solid line shows the fo contribution with Ct° = 9.3 GcV -1 . (a) shows 
the result in the entire energy region, and (b) shows that in m^o^o > 850 MeV. Experimental 
data shown by o are from Ref. |2Zj, and those by • are from Ref. |28| . 

3 Charged K loop contribution 

Now let us explore the i^-loop contributions to the radiative <j) decays. The relevant Feynman 
diagrams are shown in Fig. El The diagrams (c) and (d) each give the same result while (a) 
and (b) are required by gauge invariance. Notice from Eq. ()2.3|) that the direct photon-two 
pseudoscalar vertex vanishes in this model when k = 2 is adopted, as we are doing here. * 3 
Thus the two pseudoscalars first couple to p, uj and 4> which then transform to a photon as 
shown in Figs, (c) and (d). The strong vector- two pseudoscalar interaction vertices may be 
read from the fourth term of Eq. ()2.2|) while the scalar-two pseudoscalar interaction vertices 
are derived from the A, B, C and D terms of this equation (and explicitly given in Eq. ()2.9j0 . 

Note again that the Lagrangian density of Eq. ()2.2|) treats all of the pseudoscalars, scalars 
and vectors in a consistent chiral invariant manner. It can be modified to include gauge 
invariant photon interactions by making the replacements: 



# 3 In the present analysis, we just use k = 2 for simplicity in calculation so that two kaons couple to gamma 
only through vector meson intermediate lines keeping vector meson dominance. 
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Figure 6: Feynman diagrams for the charged K loop contributions to <fi(p, e v ) — > ao(Q)+j(k, e). 
The solid line denotes the ao meson, the wavy line the photon, the double solid line the vector 
mesons (p, u, 4>) and the dashed line the K meson, p, k and Q are the momenta of the (j> 
meson, the photon and the do (980), respectively, and I is the loop momentum. 
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where A and Q were defined after Eq. (|2.3|) . Under an infinitesimal electromagnetic gauge 
transformation with 5A^ = d fl X(x), p^ and — gp^ in Eq. (|3.1|) do not contain any terms 
proportional to d^\(x). When substituted into Eq. (|2.2|) . the above replacements yield, in 
addition to Eq. (j2.3J) the four field photon interaction terms in the Lagrangian density: 

^A^K+K- + l e J^L A ,al(K + d,K~ - K~d,K + ) + ■ • ■ , (3.2) 

V2 

where M is the 0- meson field and Oq is the neutral ao(980) scalar meson field. * 4 Now 
it is straightforward to obtain the K loop amplitudes (with the assumption k = 2) for 
0(p,e v ) — ao(Q)+7(M: 

dH [l-(Q-l))[e-e v ] 



S a = hJ 
S h = 



z(2vr) 4 [l 2 + m 2 K ][(Q-l) 2 +w? K y 
h r dH [(21 -p) -e v }[(2l-p) -e] 



2 7 z(2vr) 4 [/ 2 + m2 f ][(p-/) 2 + m 



A' J 



a c h r dH [l.(Q-l)][(2l + k).e][(2l + k-Q).e v } 
c d ~~ 2 J i(2nY [I 2 + m\\ [(Q - I) 2 + m\\ [(k + I) 2 + m 2 K ] ' 1 ' ' 



where 



h = 6T7 y a ff « V2eg^ aKK , (3.4) 
V2c/F 2 

and the KSRF relation was used in the last step. Note that the quantity defined in Eq. (|2.8|) . 
g 2 = — Q 2 . To get the amplitude for the decay — > / 7 we should replace j a KK by ^jkk in 
Eq. (jS3D. 

The next step is to regulate the divergences which occur in these amplitudes. We employ 
the dimensional regularization scheme and thus continue from 4 to d space-time dimensions 
according to the formula: 

d d l 1 r(n - d/2) 



i(27r) d (I 2 + s) n {A<K) d l 2 T{n)s n - d / 2 ' 



(3.5) 



where n is an integer while s is arbitrary. The physical amplitudes will emerge in the limit 
when e = 4 — d ^ 0. It is convenient to define: 



# 4 The terms such as e^edefN^gp^ — v^)^(gp^ — v^)) can be added into the Lagrangian (|2.2|) . Although 
they do not contribute at tree level to the radiative decays studied in the present analysis, they generate the 
vertex of type SVjPP, which gives the quantum correction to SV"f vertex. Since this quantum correction 
does not depend on the external momenta, its contribution is absorbed into the redefiniton of the effective 
SVj coupling Cy. 



14 



1 2 

- = - -7 + ln(47r), (3.6) 

where 7 ~ 0.577 is the Euler-Mascheroni constant. 

For S a we use the identity l-Q-l 2 = -l/2[(/-Q) 2 + m 2 K ] - l/2(/ 2 + m\) + [m\ + Q 2 /2) 
to write: 

S a = ih5^el(p)e v (k)[-A (m 2 K ) + \{2m\ + Q 2 )B (Q% (3.7) 



where 



d d l 



i{2ix) d l 2 + m 2 ' 



, d d l 



Bo{Q ) I i(2<n) d [l 2 + m 2 K }[{l-Q) 2 + m 2 ' ' ^ 



For Sb we define: 



MPj 7 z(27r) d [/ 2 + m 2 ,] [(p - I) 2 + m\\ ' 1 ' 

Finally, for the triangle diagrams we similarly rearrange the numerator to get: 



S c = S d = -i^el{p)e v {k)[~B^(k) + Q^B u (-k) + ^Q,KB (k 2 ) 



^£Up) - B,{p)Q v + l -p,Q u B Q {p 2 ) + 1(277-4 + Q 2 )^fe *0L (3-10) 



wherein B Q (p 2 ) and B^ u {p) have been already defined while, 



B ^ P) I i(27ry{l 2 + m 2 K ]{(l-p) 2 + m 2 K y (3 ' H) 



and 



Note that k 2 = since it corresponds to a physical photon momentum. 

Using Feynman's trick for combining denominators and Eqs. ()3.5j) and ()3.6)1 we evaluate 
the integral B (p 2 ) near d = 4: 

5o(p 2 ) = 7 ^[i-Fo(p 2 )], 

F (p 2 )=[ dx\n[m 2 K -x(x-l)p 2 ]. (3.13) 

JO 
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We also find 

5» = \p»B,{p 2 ). (3.14) 

The presence of a pole at d — 4 indicated by the term 1/e corresponds, of course, to a 
logarthmic divergence in the cutoff regularization method. 

Using the integrals defined above we can compactly write the total amplitude as: 

S = S a + Sb + S c + Sd 

= i^el(p)e v (k)[5^(-2A (m 2 K ) + (2m 2 K + Q 2 )B (q 2 )) + B^k) 

- (2m 2 K + Q 2 )X^(p,k)]. (3.15) 

Notice, in particular, that the contribution of Sb has cancelled out against a piece of the 
triangle diagrams. 

The evaluation of an integral of the form B^ v (p) is a little more complicated. We use 
covariance (in (^-dimensions) to relate it to the other integrals as: 

+ pM^({ 1 + ^r^ )Ao(m ^ - m ^ (1 + T^d )B ^ p2) 
-^Y^d Bo{p2)) - Bo{p2)] - (3 ' 16) 

We see that B^ contains the integral v4o(m| : ) which is noted from Eq. ()3.8j) to involve a 
quadratic divergence in the cut-off regularization scheme. In the dimensional regularization 
approach this corresponds to a pole at d — 2, as may be seen from Eq. ()3.5|) . It is 
necessary to check that this divergence cancels out in the total amplitude. This may be done 
by using Eq. (|3.1fij) to get, near d = 2: 

B„ u (k)\ k 2 =0 = 2A {m 2 K )5^ + ■■■ , (3.17) 

where the three dots indicate terms not containing A (m|-). Substituting this into Eq. ()3.15j) 
(considered at d=2) shows that all dependence on A (m|-) at d = 2 is cancelled, as desired. 
We interpret this as the cancellation of the quadratic divergences in the individual diagrams. 

For the physical case we must consider, of course, the amplitude evaluated near d = 4. 
The integral A (m 2 K ) is, near d = 4: 



A (ml 



K) 



1 , . 2 

- + 1 — m(m K 



(An, 

Using Eq. ()3.16j) we find for B liv (k) near d = 4 and k 2 = 0: 



(3.18) 
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B^{k) = 2A {m 2 K )8 fMV + 



1 w 2, 2 
- - \n{m K ) + - 



(3.19) 



3(4tt) 2 m " 

wherein the first term was separated for convenience. Note that the k^ky term will not 
contribute to the physical amplitude because it gets multiplied by the photon polarization 
vector e u (k). Now substituting Eq. ()3.19j) into the total amplitude, Eq. ()3.15|) . shows that 
its effect is simply to cancel the — 2Ao(m 2 < ) term. 

To evaluate the remaining, X^lp, k) term we first use covariance to express it as: 



X^(p, k) = 5^X 1 + Pf,p u X 2 + k^Xz + p^KX^ + k^ v X h 



(3.20) 



where each of the Xi depends on p 2 and p-k. The Xi may be determined by calculating X^, 
kfj.X^, PyX^y and k^p v X^ u both from Eq. ()3.20|) and from Eq. (|3.12j) . This leads to the 
relations (remembering k 2 = 0): 



X 2 
B (Q 2 ) 


x 5 



0. 

X 1 + k-pX 5 , 
k-pX 3 +p 2 X 4 , 



1 



p- 



+ 



(p-k) 2L d~2 
Am 2 K k-p 



[B (Q 2 ) - B (p 2 )) + (p-k 



Q 2 -P 2 
d-2 



)B (Q 2 



d-2 



C(p 2 ,k-p)}, 



(3.21) 



where the finite integral C(p 2 , k-p) is given by: 

d 4 l 



C(p 2 , k-p) = J - 



(3.22) 



i(2vr) 4 [P + m\\ [{I - Q) 2 + m 2 K ) [(I + k) 2 + m\\ ' 

Actually, only the coefficients X\ and X 5 remain after X^ u is multiplied by the polarization 
vectors of the photon and 4> meson; furthermore these two coefficients are related as above. 
Substituting back into the total amplitude, Eq. (|3.15|) and making use of the cancellation 
between the Ao(m 2 K ) and B flv (k) terms discussed before, yields: 



S 



V - {B (Q 2 ) - B (p 2 )} 



+ 2m 2 K (p-k)C(p 2 ,p-k) 



p-k 



(3.23) 



Note that the last term arises from the 1/e term in B (Q 2 ) multiplying the leading e term of 
its factor. From Eq. (J3.13)) we see that the logarithmic divergences cancel out of the difference 
(B (Q 2 ) — B (p 2 )). Thus the final amplitude is completely finite; both the logarithmic and 
quadratic divergences have been seen to cancel using regularized expressions for everything. 



17 



The quadratic divergences arose in the first place because of the derivative-type interactions 
required by use of the non linear sigma model terms to describe the pseudoscalar meson 
interactions. In addition, the starting Lagrangian treated the vector and scalar mesons in 
the same chiral invariant framework. 

Evaluation of the finite integrals in Eq. ()3.23J) yields the final expression for the Feynman 
amplitude, iS: 



where, 



iS 



h 1 



<(pM*)(- 



+ 



p 



2 (An) 2 " v /v p-k (p-k) 

1 1 l + a(p 2 



k -^)(2m 2 K + Q 2 ){-p-k 



m\ 



a(Q 2 )[\n 



+ 



m 



K 



l + a{Q 2 ) 

i - °(Q 2 ] 



— m\ 



[In 



l + a(p 2 ) 
1 - a(p 2 ) 



n 1 + g W 2 ) ■ v 



(3.24) 



<Q 2 ) 



\ 



i + 



Am 2 K 



\ 



1 + 



Am? K 



Q 2 



(3.25) 



Note that Eq. ()3.24|) holds only in the kinematical range where — Q 2 = q 2 > 4m 2 K ; the 
positive quantity, q was also defined in Eq. (|2.8[) . Furthermore note that p 2 = —m 2 ^. In the 
kinematical range where — Q 2 = q 2 < Am 2 K) one should replace 

l + a(Q 2 ) 



[In 



«7T 



l-a(Q 2 ) 
a(Q 2 ) - *a(Q 2 ) 



— 2itan 1 



1 

HQ 



2V 



Am 2 K 



Q 2 



(3.26) 



in Eq. (|3~2~4"j) above. 



4 Comparing the loop with experiment 

The expression in Eq. (13.24)1 describes the decay <p a ol- To get the Feynman amplitude for 
(f> — > 7i°ru, we should multiply Eq. ()3.24j) by the factor (qi ■ q2)D ao (q 2 )'j ao7rri , where D ao (q 2 ) 
was defined in Eq. (12. 7j) . This assumes a simple form for the do propagator, which can 
only be an approximation. However our main concern here is an initial exploration of the 
resonance region in the present framework so it seems reasonable for now. We write the 
resulting Feynman amplitude as 
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,5(0 _> n o m) = e(gi . q2)x (, v) [ (p . k){( y . t) _ {p . t){k . e V)j > (41) 

which thereby defines Note that the sum Xfw) + Y^ v \ where Y^^ is defined in 

Eq. ()2.6|) . would correspond to a model containing both the K loop contribution to the 
resonant amplitude as well as a point vertex contribution to the resonant amplitude. For 
now we focus on the i^-loop contribution. The decay spectrum shape, dT/dq is then obtained 
by replacing Y^) [ n Eq. (EHUl) by Xj™\ Conventionally, one uses instead, 

dB((f> -> 7rV) _ 1 dT{<P -> TT ^) 

<ig r(0) <ig 

where = 4.26 Me V. 

In section 12.11 we observed that, even though the use of derivative type SPP coupling 
helped somewhat, the tree interaction involving the ao(980) resonance was unable to explain 
the shape of the peak at large q in the experimental data for dB(<p — > ^rj^/dq. Now 
we will look at the result of using the K-\oop amplitude for this purpose. Taking [20] 
m a = 985 MeV and T a = 50-100 MeV, the only quantity which is not well known is the 
product of the scalar meson coupling constants "faKKla-nri- 111 Fig- El h is shown that a choice 
JaKKlawr] = 125 GeV~ 2 can nicely explain the shape of the experimental data in the region 
of q near the ao(980) resonance. 
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Figure 7: Predicted dB((f> — > ir°r]"f)/dq in the region of the ao(980) resonance with ^aKKla-K-q = 
125 GeV~ 2 and T a = 0.1 GeV. The vertical scale has units 10 _7 MeV _1 . Experimental data 
are as in Fig. 

For q below the resonance region, the K loop contribution in the present model falls off 
rapidly, as one might reasonably expect with derivative coupling, and lies lower than the 
data points. In addition to the nonresonant background [3] which is usually included to 
explain this region, there might be some tree level resonance production which was observed 
in Fig. 121 to peak around 950 MeV. 
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The main feature of the data is that there is a very rapid falloff with apparent discontinu- 
ity of the slope, when q reaches the KK threshold. This is a clear signal for the importance 
of the K loop contribution. One may see this feature by referring to Fig. |HJ for which the 
ao(980) mass has been artificially lowered to 970 MeV. Comparing with the previous figure 
shows that the sharp fall-off is exactly the same in both cases, clearly unaffected by the dif- 
ference of assumed resonance masses in the two cases. The difference in masses, on the other 
hand, shows up as a difference in position of the peaks. It should be remarked that the peak 
position is also affected by the decreasing phase space with increasing q. This characteristic 
feature of the iT-loop contribution was first illustrated by Achasov [30 by considering the 
behavior of the result with lowered values of the i^-meson mass. 
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Figure 8: Predicted dB((f> — > ir°r)"f)/dq in the region of the ao(980) resonance but where the 
ao mass was artificially lowered to 970 MeV. Here jaKK^a-nrj = 115 GeV -2 and T a = 0.1 GeV. 
Experimental data are as in Fig. [21 

Next, let us check the dependence of the prediction on the width of a (980). Figure El 
shows that the predicted dB(<p — > 7r°^7) / dq in the region of the a (980) resonance with 
laKKla-wri = 95GeV~ 2 and Y a = 0.05 GeV. Comparing this prediction with that given in 
Fig. [3 we see that the smaller ao width gives a sharper peak, and that a smaller value 
of laKKla-Kri can also reproduce the experimental data at the peak position. For further 
decreasing the value of •jaKKlanri the inclusion of the K-loop correction into the propagator 
of ao(980) may be important as pointed out in Ref. [30J. 

The i^-loop contribution to the branching distribution, dB{<p — > 7i TT 0/ ~f)/dq may be 
similarly evaluated and compared to experiment. There is similarly a problem for the tree 
level resonance model to reproduce this experimental shape in the high q region. The K loop 
amplitude <fi — > fo (980)7 is given by Eq. (|3.24|) wherein the overall factor h is now obtained 
by replacing j a KK in Eq. (|3.4|) by •jfKK- To get the Feynman amplitude for — > 7r°7r°7, we 
should multiply Eq. (J3.24|) by the factor V2(q% ■ q2)Df Q (q 2 ) / yf 07r7T , where Df Q (q 2 ) was defined 
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Figure 9: Predicted dB(tfi — > ir°r]"f)/dq in the region of the ao(980) resonance with jaKKjaTrr) = 
95 GeV -2 and T a — 0.05 GeV. Experimental data are as in Fig. |2 

in Eq. (J2.19j) . This defines X^ as in Eq. (|4.1j) . The spectrum shape is determined by using 
Eq. (J2~T7J) with X ( / o n) replacing . 

Taking [20] m f = 980 ± lOMeV and T f = 40-100 MeV, the only quantity which is not 
well known is the product of the scalar meson coupling constants "ifRKlf-mr- In Fig. [TUJ it is 
shown that a choice jaKKlf-Kn = 86 GeV -2 can nicely explain the shape of the experimental 
data in the region of q near the /o(980) resonance. 
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Figure 10: Predicted dB((f> — > ir°Tr°j)/dq in the region of the /o(980) resonance with 
7/KK7/7T7T = 86 GeV" 2 and T f = 0.1 GeV. The vertical scale has units 10" 8 MeV _1 . 

As in the case of the <fi —>■ 7r ?77 process, the i^-loop description of dB/dq only explains 
the upper q region near the scalar resonance. To cover the lower q region some non resonant 
background [3] is required. Possibly a tree level resonant background, corresponding to using 
Xj*^ + in Eq. (|2.17|) . would also be appropriate. 

In both cases considered in this section, it is also desirable to include the effects of 
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using multichannel scalar meson propagators |3T| for a better approximation to the detailed 
dynamics. A very recent treatment of the — > 7r°7r°7 process in this framework is given in 

5 Summary and discussion 

Historically, the study of elementary particle spectroscopy has been built around the orga- 
nization of these particles into SU(3) flavor multiplets and the consequent predicted (broken 
symmetry) mass formulas and interaction vertices. The still mysterious scalars can be ex- 
pected to yield up some of their secrets by this type of analysis. Indeed some recent analyses 
have already been carried out [TTJ |33J EU EH]- The most dramatic feature is that the 
light scalars appear to exhibit, as originally suggested by Jaffe [37], a reverse mass ordering 
compared to the other meson multiplets. 

In Ref. jU] an attempt was made to extend the SU(3) analysis to relate all the decays of 
the types S— > 77, V-^S7 and S— 3V7 to each other by using a simple effective VVS point- 
like interaction, together with vector meson dominance. The analogous assumption 7J of a 
point-like VVP structure was very successful [HH] i n phenomenologically correlating P— > 77, 
V^P7 and P— 3V7 decays. Such an approach was the original motivation which led to the 
present investigation. In section 2 we compared the spectrum shape of the decays </> — > r/7r°, 
measuring the effects of an intermediate ao(980) resonance and <fi — > 7r°7r°, measuring the 
effects of an intermediate /o(980) resonance, in the point like VVS model with the corre- 
sponding experimental observations. It was found that the resonant peaks in the model were 
pushed lower due to decreasing phase space. This contrasted with experiment which does 
not indicate this effect. On the other hand, if one were to use a tree model of this type with 
non derivative SPP type couplings, the resonant peaks were seen to get completely washed 
out. This would appear to be an advantage for the derivative coupling, which is dictated by 
chiral symmetry in the present framework. Nevertheless, since even with derivative coupling 
the spectrum shape is not very well fitted, there must be another mechanism at work. 

Now, it has been emphasized [30J that the f^-loop model for the <fi radiative decays 
constitutes a special mechanism which does give a characteristic spectrum shape in agreement 
with experiment. This is readily understandable since the 0(1020) meson is just a little bit 
heavier than the two K mesons which comprise its main decay product. We thus studied the 
.fT-loop diagrams using the chiral Lagrangian of pseudoscalars, vectors and scalars given in 
Eqs. (|2.2|) and (|2.4j) with the relevant photon terms introduced by the substitutions shown in 
Eq. (|3.1|) . Most of the calculations of this process have not started from a chiral symmetric 
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Lagrangian and have thus used non- derivative type SPP type interaction vertices. The use 
of derivative coupling introduces an extra complication in that there is an new diagram, 
shown as (b) in Fig. |U1 In addition, individual diagrams now contain quadratic as well 
as logarithmic divergences. It is known that these divergences are forced to cancel from 
gauge invariance. However we have used the dimensional regularization scheme and shown 
explicitly in section |H] that both the log and quadratic divergences cancel in the regularized 
expressions. This may be of some interest in dealing with processes of the present type. 

In section |H we observed that the shape of the ao(980) and / (980) resonance regions 
in the radiative decays could be explained by the corresponding K-loop amplitudes. Fur- 
thermore, it was evident that the characteristic sharp drop in the amplitude at large q 2 
was associated with the the KK threshold rather than with the falloff of the resonance 
away from its peak. For this work, we used the coupling constant products ■jaKK'Janr] and 
IfKKlf-K-K respectively as fitting parameters for the radiative decay spectra into tt i] and 
7T°7r . Elsewhere, we plan to study more precisely the values of these coupling constants 
obtained by comparing with experiment, chiral models of meson meson scattering in which 
the same interactions are used. We also will study how the point-like diagrams with resonant 
contributions discussed in section El can be used in conjunction with the i^-loop diagrams 
to improve the fit to the resonant region. This will presumably become even more interest- 
ing when more data points become available. Another point of interest concerns the extent 
to which the various SPP coupling constants can be correlated assuming a single nonet of 
scalars. This arises because there is some evidence that two scalar nonets (one presum- 
ably made from 4 quarks and the other from two quarks) mix to make up the physical scalar 
states. A recent exploration of the effect of such a mixture on 0(1020) radiative decays has 
been given in Ref. jlHI- Still another correction to the simple picture employed here would 
be to use more realistic resonance propagators by including pseudoscalar loops pH] . 

Of course, in order to make a careful comparison with experiment one should include 
non-resonant contributions which are expected to dominate for small q. These will include 
the emission of a pion with a virtual p which subsequently decays into 7T7 (and similar 
diagrams leading to a 7r°?77 final state) as discussed in Ref. j^j. There will also be non 
resonant contributions from the K-loop diagrams. A variety of interference mechanisms to 
explain the full spectrum are discussed in Ref. [41j. It should be noted that the "background" 
contributions may very well have a non-trivial effect also in the resonance region itself. 

Using the results obtained here and taking into account the features just discussed, we 
will continue to study the radiative decays with the expectation that it may contribute to 
the understanding of the puzzling scalar mesons and ultimately to low energy QCD. 

DEDICATION: We are pleased to dedicate this paper to Rafael Sorkin in connection 
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with the world wide web celebration of his sixtieth birthday. We have benefited a great deal 
from our interactions with Rafael. We wish him good health and continued success in his 
endeavor to understand the deepest mysteries of space-time structure. 
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